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We study the finite volume/temperature correlation functions of the (l+l)-dimensional SU(W) principal chi¬ 
ral sigma model in the planar limit. The exact S-matrix of the sigma model is known to simplify drastically 
at large N, and this leads to trivial thermodynamic Bethe ansatz (TBA) equations. The partition function, if 
derived using the TBA, can be shown to be that of free particles. We show that the correlation functions and 
expectation values of operators at finite volume/temperature are not those of the free theory, and that the TBA 
does not give enough information to calculate them. Our analysis is done using the Leclair-Mussardo formula 
for finite-volume correlators, and knowledge of the exact infinite-volume form factors. We present analytical re¬ 
sults for the one-point function of the energy-momentum tensor, and the two-point function of the renormalized 
field operator. The results for the energy-momentum tensor can be used to define a nontrivial partition function. 
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I. INTRODUCTION 


One of the main goals of statistical physics is to calculate expectation values of observables in a system at finite temperature. 
The partition function usually contains enough information to find some of these expectation values, which can be computed 
by taking different partial derivatives of it. The most common tool used to derive the finite-temperature partition function of a 
two-dimensional integrable quantum field theory is the thermodynamic Bethe ansatz (TBA). In this paper we propose that this 
tool does not work in general for a matrix-valued quantum field theory. The expectation values of operators include additional 
information not contained in the TBA partition function. 

To demonstrate this proposal, we examine the (iH-l)-dimensional principal chiral sigma model (PCSM). The PCSM has the 
action 


5pcsm= j cPx^Tvd^.U\x)d^U{x), (I.l) 

J ^9o 

where U(x) G S\J{N). This model has been shown to be integrable, and its exact S-matrix is known iH . The action ( II. Il l has an 
S\J{N) X S\J{N) global symmetry given by U{x) -A VpU{x)Vii, with € SU(N). The PCSM is asymptotically free and 
has a mass gap, which we call m. In our analysis, we simply assume the existence of a mass gap. A mechanism explaining how 
this mass is dynamically generated was proposed in ||2l. 

We are interested particularly in’t Hooft’s large-W limit of the PCSM. In this limit, the S-matrix greatly simplifies. We later 
show that this means that the TBA partition function at large N is that of a free theory. However, we show that the expectation 
values of operators are not trivial. 

The inefficiency of this partition function is due to the matrix structure of the theory. It has been shown that the form factors 
of operators are not trivial at large N iMil. As we will see, the computation of thermal expectation values can be done by 
summing over form factors. Since these form factors are not trivial (despite the trivial S-matrix), the thermal expectation values 
are not trivial either (despite the trivial TBA). 

In the rest of this paper we will use interchangeably the terms finite volume, and finite temperature. This is because in l-i-l 
dimensions, these two are equivalent up to a Wick rotation. 

In the next section we review some of the exact results that are known for the PCSM. In Section III we show how the form 
factors have been used before to calculate infinite-volume correlation functions of operators, which agree with the asymptotic 
freedom of the PCSM. We discuss the results of Ref. |@], where the infinite-volume two-point function of the renormalized-field 
operator was computed. 

In Section IV, we discuss the application of the TBA to the PCSM. We find that the TBA yields the partition function of a free 
field at large N. 

In Section V we compute the one-point function (vacuum expectation value) of the trace of the energy-momentum tensor 
operator. This correlation function is computed using the so-called Leclair-Mussardo (LM) formula H. We observe that this 
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one-point function does not agree with what is expected from the trivial TBA. We then show how to define a nontrivial partition 
function from our result for the energy-momentum tensor. 

In Section VI, we compute two-point function of the renormalized-field operator in a finite volume. This is a finite-volume 
version of the result of iQ]. While the LM formula for one-point functions is generally believed to be accurate, the validity of the 
LM formula for two-point functions has been disputed H, d. We argue that the objections from and do not affect the 
PCSM at large N, and that the LM formula might be valid in our case (though we have no proof that this is the correct two-point 
function). The very large and very small volume limits of this two-point function are examined in detail in Section VII. 


II. FORM FACTORS OF THE PRINCIPAL CHIRAL SIGMA MODEL 


In the following two sections we show a brief review of previous results on exact form factors and correlation functions of the 
PCSM. 

The main tool that has been used in previous works is the form factor bootstrap program for integrable field theories ifl^ . The 
integrability of the PCSM implies that all scattering events are completely elastic and factorizable into a product of two-particle 
S-matrices. 

All the qualities of an elementary excitation are specified by stating its rapidity 6, related to its energy and momentum by 
E = m cosh 9,p = m sinh 9, its left and right SU(A^) color indices a,b = 1,... ,N, respectively, and by stating if the excitation 
is a particle or an antiparticle. We can write, for example, a one-particle incoming state and a one-antiparticle incoming state as 


\P^ 9^ U, b'^in: \^: ^)in5 


respectively. 

The particle-antiparticle S-matrix, defined by 


out(^, 0'i, di, Ci; P, 0', C2, dM, 01, 5i, ai; P, 02, 02 , 62 )i„ = " d'l) 4^<5(02 - 0^), 


is known to be lH 


o(n\d2C2\cidi 

^ ^aibi-,b2CL2 


Q{0) 


AC1AC2__A AC1C2 

Af(7ri - 0) 




27ri 


A^(7ri — 0) 


Jbib2 




where 


Q{0) = 


sinh 


2 


sinh 


2 


7ri 

N 


El 

N 


( r[j(7ri — 9)I2tt + l]r[— i(7ri — 9)(2tt — 1/N] 1 ^ 
\ r[i(7ri — 0) /2tt -f 1 — 1 /A^]r[—i(7ri — 0)/27r] J ’ 


(II. 1) 


( 11 . 2 ) 


and 0 = 01 — 02. The particle-particle and antiparticle-antiparticle S-matrices can be found using crossing symmetry. An 
incoming particle (antiparticle) can be turned into an outgoing antiparticle (particle), by shifting its rapidity by 0 —0 — ttL For 
general N, there exist r-particle bound states, with mass 


Ulr 


sin(f) 

m- 

sin(f) 


r = 1,...,A^- 1. 


For the rest of this paper we will work exclusively in ‘tHooft’s large-W limit. That is, we take N — 00 , while keeping 
m fixed.This limit simplifies the problem in many ways. First of all, there are no bound states at large N, since the binding 
energy vanishes. Also the S-matrix is greatly simplified in this limit, as Q{9) = 1 + 0 (l/iV^). As was pointed out in ifTl]] . 
the thermodynamic Bethe ansatz equations in ‘t Hooft’s large-A^ limit are essentially those of a free theory (the authors of this 
reference later investigate a different large-A^ limit with nontrivial Bethe equations, which we do not discuss further here). 

At large N, two excitations interact nontrivially only if they have color indices contracted with each other. This is easily seen 
from Eq. (HLB. The non symmetric terms in the S-matrix, proportional to or Sbib 2 d‘^^'^^^ vanish at large N, unless 

one sums over the colors of one of these delta functions. A particle has a left and a right color index, so it can interact nontrivially 
with at most two other excitations. 

The form factors (matrix elements of local operators) of the renormalized field, $ (a;), have been found in the large-limit, 
in Ref. JH]. This field is defined in terms of the bare field, U, by 

(0|Tr$(x)$(0)t|0) = Z[ 5 o(A), A]-i(0|TrP(x)P(0)t|0), 
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where Z[( 7 o(A), A] is a renormalization constant, A is the Euclidean momentum cutoff, and ( 70 (A) is the coupling constant, 
which runs such that the mass gap, m, is independent of the cutoff. We write here form factors with excitations only in the 
incoming state, as outgoing particles can be obtained using crossing symmetry. Because of the SU(A^) x SU(7V) symmetry of 
the PCSM, only form factors with M particles and M — 1 antiparticles are non-vanishing, where M is a positive integer. The 
form factors can be parametrized as 


(0|<l>(0)b(,a(, \A, 01, &1, oi;...; A, bM', ■ ■ ■ ',P,&2M-i,a2M-i,b2M-i) 

= iV-^+1/2 ■ ■ ■ , d2M-l) 


(11.3) 


where cr is a permutation that takes the set of numbers 0,1, 2,..., M — 1 to (t( 0), cr(l),..., a{M — 1), and r takes the numbers 
0,1, 2,..., M — 1 to r(0), t(1), ..., r(M — 1), and we sum over all the possible permutations in the set Sm- The main result 
of Ref. is (at large TV): 




. . . , 02M-i) = < n3 = i [Sj-8„ij)+M+Tri\[6j-e-rU)+M+T^iV 

1 0 , otherwise 


o-(j) 7 ^ T{j), for all j 


(11.4) 


A crucial tool for being able to find these form factors was the simplicity of the S-matrix at large N. The scattering of any two 
incoming excitations in ( III. 3b is trivial except for the pairs of permutations cr, r where one or both of their color indices are 
contracted. 

We are also interested in the form factors of the energy-momentum tensor. These have been found in Ref. a. By SU(N) 
color symmetry, non vanishing form factors have the same number of particles and antiparticles. These are 


(0 I T^jy( 0 )|A, 01, 61 , tti;...; A, 9m, bM, um', P, 9M+i,aM+i,bM+i', ■ • ■; ^ 2 M, o,2m, b2M) 

— [{pi + • • • + P 2 m)ii(P 1 + • • • + P2m)v ~ Pfj,v{pi + • • • + P 2 m)^] 

^ MM 

^ _/VM-l E ( 01 , ... , 02 m) n ^Chj 0 'cr(j) +Af x(fc) + M ; 

k—i 

where a and t are the permutations that take the numbers 1,..., M to cr(l),..., t{M) and t(1), ..., t{M), respectively. At 
large N: 


( —27r^)(4'7r)^ 


1 0, otherwise 


, for cr(j) ^ r(j), for all j. 


(11.5) 


We will be interested in the trace of the energy-momentum tensor operator, 0 = T)). 


III. CORRELATION FUNCTION AND ASYMPTOTIC FREEDOM 

Using the exact form factors. Eg. (III. 31) and (III. 4b . an expression for the infinite-volume two-point correlation function of the 
renormalized field was written in Ref. (j]. This correlation is found by summing over all the intermediate states: 

^ ^ E Ee“'""(0|‘J’(0)boaol^)('I'l[‘i>(0)boao]10), (III.l) 


ao,6o 


ao,&o ^ 


where 'k is any state with particles and antiparticles, andp^ is the sum of the momenta of the excitations of the state T*. 
By directly introducing the exact form factors into (llll.lb . one finds 


W(x) = ^ 


M 


1 1 
(M- 1)! M! 


/ 2M-1 




n ^ ) El^‘"^(^i’---’^2M-i)pexp j Ta:- ^ Pj 


2M-1 


-o 


i=i 


A drastic simplification comes from realizing that for a given value of M, all the pairs of permutations cr, r give the same 
contribution to the correlation function. The final result of Ref. il is (ignoring 0{1IN) terms) 


2M-1 


W(2^) = ^E/ d9i... (i 02 i+i exp ia: • ^ p J PJ —— 


1=0 


i=i 


Ll - dj+i)"^ + 7 r 2 ■ 


(111.2) 
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Recently, the short-distance behavior (x 0) of the function (IIII.21 I has been examined |@]. This was done in Euclidean space 
by looking at =0 and = iR, such that expja; • Pj —>■ exp —mi? cosh 0^. The strategy is to realize that for small R, the 
function exp —mR cosh 9j looks like a plateau, where it is approximately 1 for —L < e, < L, and zero everywhere else, where 
L = In This technique was first used to study the short-distance behavior of the Ising model ifl^ . 

For short distances the function (I1II.21 I then becomes 


W(ii?,0) 



1 

0j+i)2 +Tr‘^' 


(III.3) 


The function (IIII.3I I was studied in i^, and it was shown that it diverges in a way that is consistent with what is expected from 
asymptotic freedom. 

There is an alternate (and equivalent) way of examining the short distance behavior of the correlation function. The function 
W(a:) diverges at the point x = 0 when one performs the integrals over the rapidities. One solution is to simply introduce a 
cutoff in the rapidities. A, “by hand”. One then finds 


>V^(0) 



1 

+ + TT^ 


(111.4) 


The function (IIII.41 I is exactly the same as (IIII.31 I. except we have replaced L by A. It is convenient to introduce the variables 
Uj = 9j /A, such that 


OO «i 


21 


= I du2i+iY[ 


1^0 


1-1 


1 Kiuj - Uj+i)^ + 7rVA2] ■ 


We define the function 


T{uj,Uk) = 


1 


A[(Ui - ■Ufc)2 -I- 7r2/A2] ’ 

and an operator T and vector space {|it)}, such that 

{ui\f\uk) =T{ui,Uk), {ui\uk) = S{ui - Uk), 1= / duj\uj){uj\. 


(III.5) 


( 111 . 6 ) 


In terms of the operator f, the correlation function can be written as 

W^(0) = E I du'du{u'\^^\u). (III.7) 

The technique used in |@] was to realize that the operator T can be written approximately in terms of the fractional Laplacian 
operator, —d?-/du^. It is shown in |0] that one can write 

f = 

where H{\) is some operator that satisfies H{\) = + 0{1/X). The fractional Laplacian satisfies the eigenvalue equation 

= ctnPniu), where n = 1, 2,... and 0 < ai < 02 < ..., with (^„(±1) = 0. 

The correlation function ( IIII.SI l can be written, for large A, as 


= IT. 


n—1 


duipn{u) 


'-1 


1 


1 _ g-277a„/A-|-0(l/A2) 


87r2 


E 

n—1 


duifniu) 


(III.8) 


The correlation in (IIII.8I 1 is proportional to A^. This rapidity cutoff is related to a standard Euclidean momentum cutoff. A, by 

sinh^(A)-I-cosh^(A) = A^, 


so that 


A = sinh 


-1 






This dependence of the correlation function on the logarithm squared of the momentum cutoff is a confirmation of the asymptotic 
freedom of the model, and is predicted by perturbation theory ifT^ . 
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IV. THE PCSM AT FINITE VOLUME 

We make the direction finite by imposing periodic boundary conditions + U) on all wave functions, where 

V is the one-dimensional volume. Placing the system in a finite volume discretizes the energy spectrum. The quantization 
condition depends on the exact S-matrix, and is found using the Bethe ansatz. This is, for an n-excitation state iGl 

n 

e^p^^Y[S{9,-9k)=±l, j = l,2,...,n, (IV.l) 

where we have suppressed all the color indices in the S-matrix for simplicity. The selection rules are ± = -f, for boson-like 
interactions, S'(O) = 1, and ± = — for the fermionic case, S'(O) = —1. Equivalently, one can write 


mV sinh 9j + '^A{9j -9k) = 2 ttI, (IV.2) 


where A{9) = —i\iiS{9) and I is an integer for bosonic interactions, and a half-integer in the fermionic case. Solving the 
equations (llV.2b one can find the discrete spectrum of rapidities, 9j. 

At large N, a particle (antiparticle) can interact nontrivially with at most two other antiparticles (particles). The most nontrivial 
n-excitation state one can define is a chain of alternating particles and antiparticles, where the j-th particle (antiparticle) has one 
color index contraction with the {j — l)-st antiparticle (particle) and the {j -f l)-st antiparticle (particle). Using the S-matrix, 
Eq. (111.lb . the Bethe quantization condition for the j-th particle at large N is 


mV sinh 9j — i In f ^ 


i In 


/ +7rA 

\9jj+i-7riJ 


= 2Trl, 


(IV.3) 


where we use the notation 9jk = 9j — 9k. There are only two (instead of n — 1) terms coming from the S-matrix in Eg. (II V. 3b. 

It is useful to find the spectrum in the thermodynamic limit, where V,n ^ oo, but their ratio is fixed. It can be shown ifl^ 
that in the thermodynamic limit, the quantization condition (IIV.2b becomes 


e(0) = mV cosh(0) =F 




(IV.4) 


where (p{9) = ■^A{9), and e{9) is the so-called pseudo energy. The interpretation of the pseudo energy is that the “dressed” 
energy of a particle of rapidity 9 is given by e{9) jV in the thermodynamic limit. 

The thermodynamic limit of Eq. (IIV.3b is trivial because the first term on the left-hand side dominates over the other two 
(because there are only two terms coming from the S-matrix, instead of n — 1 terms as in the usual TBA). The TBA at large TV 
is therefore trivial. The pseudo ener gy o f a particle of rapidity 9 is simply given by t{9) = mV cosh0 (because kp{9 — 9') — 0 
in ( IIV.4b ). This was noticed in Ref. iflin . where this result is used to declare (incorrectly, as we propose) the’t Hooft limit not 
physically interesting. 

Once the pseudo energy is known, the ground-state energy, Eq{V), and partition function,Z(L, V) obtained from the TBA 
arelfll: 


/ dO 

—mcosh01og(l ± 
27r 

Z{L,V) = e^v[-LEoiV)], 


where L is the size of the direction. The partition function (IIV.5b can be written equivalently as a sum over states; 


(IV.5) 


Z{L,V) 



2=1 


(IV.6) 


where the scalar products in (IIV.6b . (0i,..., 0„|0i,..., are those of a free bosonic or fermionic theory. The only effect of 
the interactions at finite volume is that the energies of the excitations are dressed. Since the pseudo energies are trivial at large 
N, the partition function derived from the TBA is that of an ideal gas. Despite this fact, it is easy to see why the expectation 
value of an operator O, is not trivial. The expectation value can be formally written as 


{O)^ = 








(IV.7) 
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The expression (IIV.7b involves a sum over the form factors of the operator. As we have seen in the previous sections, the form 
factors are not trivial, even at large N. The expectation values are then different from those of a free theory. This expression 
has singularities that need to be regularized, before it can be used explicitly. In the rest of this paper we use the regularization 
scheme proposed by Leclair and Mussardo 13]. We do not derive the Leclair-Mussardo (LM) formula here but simply quote it 
and use it. 

It is easy to understand the failure of the TBA partition function if we don’t suppress the color indices in the expression (IIV.6I) . 
The scalar products in (IIV.6I 1 involve only the symmetric part of the S-matrix, while disregarding the effect of any nontrivial 
contractions of color indices between particles. At large N, the symmetric part of the S-matrix is trivial, so all the nontrivial 
information from the S-matrix is ignored. On the other hand, the form factors in (IIV.7I) involve a sum over all the nontrivial color 
contractions. The main disparity is that the TBA partition function throws away the nontrivial color contractions, while the form 
factors do not. 

One simple modification to the partition function (IIV.6b is to use the full S-matrix elements of the form (III.lb to compute 
the scalar products, instead of using the scalar products from the free theory. The contributions from the antisymmetric part of 
the S-matrix are suppressed by higher powers of 1/N. For example, the two-excitation contribution to the partition function 
involves the scalar product 


Xaisa2^bl cb: 


ry / /r\ C 



2Tri \ ^ 

iri — 9 ) ’ 


where the TBA partition function accounts only for the leading A^'* term. The subleading terms cannot be ignored if one is 
interested in computing correlation functions. This is because only the antisymmetric part of the S-matrix gives nontrivial 
contributions to the form factors. The terms that are suppressed in the partition function become the leading terms in the 
correlation functions. 

It is important to stress out the peculiarities of the large-A^ limit that made the Bethe-ansatz calculation so simple. First, the 
large-A^ S-matrix is diagonal, as was pointed out in jjstj. When two excitations scatter, they keep their identity (color charge 
quantum numbers), as the amplitude of an identity-changing process is suppressed by a factor of 1 /N. This is clearly seen by 
examining Eq. dira . The only nontrivial scattering process is when two excitations have one or two color-index contractions. 
One can then forget about the color structure and reformulate the problem as a diagonal theory of particles that can interact with 
the S-matrices 


S{9) 


9 

’ 0 =F TTZ ’ 


/ 9 ±TTi 

=F TTZ 


2 


The second simplifying property of the large-W limit is that, as we mentioned in Section II, there are no bound states. The 
binding energy of the bound states vanishes at large N, therefore the Bethe ansatz involves only elementary particles, with the 
scattering we just described. 

The Bethe ansatz is significantly more complicated for arbitrary finite N. The particle spectrum consists of A^ — 1 bound 
states with non-diagonal scattering. The product of S-matrices in Eg. ( II V. lb (commonly referred-to as the transfer matrix) has 
a complicated color structure and needs to be diagonalized. The eigenstates of the transfer matrix can be expressed in terms 
of the physical A^ — 1 bound states, plus auxiliary “magnon” particles that carry no energy or momentum. Eurthermore, the 
Bethe equations for the auxiliary particles allow for additional magnon bound states, usually called “strings”. One can then 
express the problem as a diagonal scattering theory, but including scattering with an infinite number of auxiliary string states. 
The necessary ingredients are all the S-matrices of physical bound states and strings. A detailed derivation of the S-matrices and 
Bethe equations for the PCSM at general N is found in ifisll . 

It is necessary to understand what happens to the contribution from all these physical bound states and string states when we 
take the large-A^ limit. The contribution from the string states to the partition function at finite A^ is a result of their nontrivial 
scattering with the physical particles (The exact S-matrices can be found in Eq. (III.22) and (III.30) of ifisll ). Erom those 
expressions, it is easy to see that at large N, these string-physical particle S-matrices become 1 + 0{1/N). That is, the strings 
and physical particles stop interacting with each other at large N. There are then no contributions from auxiliary strings to the 
partition function 

The disappearance of the contribution to the ground-state energy from the physical bound states is a bit more subtle. As we 
have mentioned, the S-matrix between two elementary excitations (with no color contractions), in the large-W limit is S{9) = 1, 
yielding the partition function of a free boson. However, for any finite N, the elementary particles satisfy fermionic selection 
rules, S'(O) = —1. If one uses the finite-A^ Bethe ansatz, one has to treat the particles as fermionic. If we later take the large-A^ 
limit, we have to do so while using the fermionic rules. This process of taking the large-A^ limit after computing the ground state 
energy, instead of before, means that the S-matrix of the elementary particles will actually be given by 


S{9) = 


1 , ii9y^0 

-1, if 6» = 0 ■ 


(IV.8) 
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The TBA for a theory with the S-matrix (IIV.8I) has been examined in ifT^ . As is expected, and as is necessary for the consistency 
of our previous analysis, it is shown in that the ground state energy and partition function of the model with S-matrix (IIV.8I) 
is exactly the same as that of a free boson. 

As we have discussed before, the physical bound states dissolve into elementary particles in the planar limit. That is, the 
contribution to the ground state energy from an r-particle bound state becomes equivalent to the contribution of r elementary 
particles. However, as this bound state is dissolved, the remaining elementary particles have the same rapidity. As we discussed, 
if we take the large-limit after computing the Bethe equations, the particles satisfy fermionic rules, and therefore we must 
enforce Pauli’s exclusion principle. That is, a state with r-elementary particles with the same rapidity is not allowed, therefore 
there is no contribution to the ground state energy from the physical bound states. 

Treating the large-iV PCSM as a diagonal scattering theory from the beginning is a very useful shortcut. The more rigorous 
approach is to compute first the off-diagonal TBA and later take the large-iV limit, but this will yield the same free boson partition 
function. 


V. THE ONE-POINT FUNCTION OF T HE ENERGY-MOMENTUM TENSOR AT FINITE VOLUME 


In this section we evaluate the vacuum expectation value of the trace of the energy-momentum tensor at finite volume. We use 
the one-point function LM formula. This expectation value is interesting because it is usually easily calculated from the TBA. In 
field theories that are not matrix-valued, it has been shown that the results from the LM formula and those from the TBA agree 

iS. 

In our case, the TBA yields the expectation values of a free theory. Our position is that this is not the right value. We believe 
the value of the LM formula is the correct one, as it uses the nontrivial form factors. We believe this discrepancy is simply a 
consequence of the field being a matrix. Our approach then will be to find the expectation value of the energy-momentum tensor, 
assuming the validity of the LM formula, and then define a partition function such that it agrees with this value. 

The LM one-point function for some operator O is 

1 /* rjf) dO ^ 

(^) — ^ ^ T / ~Z ■ ■ ■ ~Z ■ • ■ 1 ; ^n)connected: (V.l) 


where f^. iOj) = 1/(1 -f (the Oj = -fl case will be relevant for the two point function). The connected form factor 

is defined as the finite part of the form factor after requiring that the rapidties of the incoming and outgoing states are equal. Any 
part of the form factors in (IV.lb that is divergent in this limit of the rapidities is discarded. This regularization is explained in 
more detail in 0 ]. 

We have used fermionic selection rules, S'(O) = —1, in defining the functions (Oj) because only the antisymmetric part 
of the S-matrix gives nontrivial contributions to the form factors. This antisymmetric part satisfies fermionic rules, while the 
symmetric part is bosonic. For a bosonic theory with S'(O) = 1, the corresponding functions in the LM formula would be 

jbosonic^^^.j = 1/(1 _ 

We need the form factors of the operator 0 = T/' with the same number of excitations in the incoming and outgoing states. 
These can be obtained from ( III.5b by crossing symmetry. 

The connected form factors of the energy-momentum tensor can be written neatly in terms of the S-matrix. We follow the 
calculation and language from Ref. 0 . For a general scalar-valued field theory with S-matrix, S{9), the connected form factors 
of the energy momentum tensor are * 

• ■ • ^n)connected — 47 rm^n! i^( 6 »i 2 )(^( 6 » 23 ) • • • i^(^'n-i.n) cosh(0i„), (V.2) 

where ip{9) = ■ 

The main difficulty when directly trying to apply Eq. ( IV.2b to our matrix-valued case is that the function <f{9ij ) is not the same 
for every pair of particles i,j. If the excitations i and j don’t have any contracted color indices, the function (p{9ij) vanishes. 
The only non-zero connected form factors are those where all the functions (p{9jj+i) are non-zero. We can build a state with 
alternating particles and antiparticles. The only color combinations that survive are those where the j-th particle has one color 
contraction with the (j — l)-st and the (j -|- l)-st antiparticles. The interaction between the j-th and the (j -|- l)-st excitations is 
given by the function 




d9 


jJ+i 


■log 


+ iri 


27r 


^L+i +' 


(V.3) 


^ This equivalence holds only for our purposes, where all the rapidites are to be integrated, as in Eq. E3. The factor of n\ accounts for the different 
permutations of the order of particle rapidities. All these permutations give the same contribution to the integral ED. 
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The non-vanishing connected form factors for our energy-momentum tensor are 


• ■ • ^a'^an ^b'^bi ■ • ■ ! P-t ^2 5 ^2 J ^2 5 ^ 3 : ^3; ^^3, . . . |0|^, 9 \^ 63, U3, O2 •> Cl 2 ; ^2; ■ • ■ ) connected 

= 47rm^n! (/?(6»i2)(/?(6'23) • • • cosh(6»i„) -f C> (i) . (V. 4 ) 


The one-point function is found by substituting the form factors (IV.41 i into the formula (IV. 11 1. Our final result is 


M! 

7V2 


= dTTTO^ ( ^ 


n / Sz-.i"*) 




n 


-27r 


eh 


i=i td+i 


cosh(0i„) 4-0 


(V.5) 


It is easy to see from Eq. (IV.4I I and Eq. (IV.5b why our results disagree with the trivial TEA. The difference between the 
connected form factors of a scalar theory (Eq. ( IV.2b ). and our matrix-valued case is that all the fundamental particles in a scalar 
theory interact with the same S-matrix. In the matrix-valued case, the S-matrix of two particles depends on how their colors are 
contracted. The Bethe equations of an n-excitation state involve the S-matrix of the j-th excitation with all other excitations. 
This is trivial in our case because excitations interact nontrivially with only two other excitations. The connected form factors 
for an n-excitation state, however, involve the S-matrix of each adjacent pair of particles j and j T 1. The n-particle state can be 
designed in such a way that all these two-particle S-matrices are nontrivial. 

Our expectation value (IV.5b can be used to define a nontrivial partition function. The expectation value of the energy- 
momentum tensor is related to the finite-volume ground state energy, Eq{V), by 


= Y^{VEo{V)]. (V.6) 

One can find the ground state energy in principle by solving the differential equation (IV.6b . The thermodynamic limit of the 
partition function is dominated by the ground state energy. We can then define the nontrivial thermodynamic limit of the 
partition function as 


Z{L,V) = 


where L is the size of the direction. 


VI. TWO-POINT CORRELATION FUNCTION OF THE RENORMALIZED FIELD AT FINITE VOLUME 


In this Section we compute the two point correlation function of the renormalized field at finite volume. For a local operator 
0{x) of an integrable theory, the LM two-point function is (again suppressing color indices) 




n=l 


dOfq 

drr 


n exp 

i=i 


(VI. 1) 


where kj is the dressed finite volume momentum of the i-th particle (which at large N is just the standard kj = m sinh 9j), and 
O is the dressed vacuum energy at finite volume. The first term in the right-hand side of (IVLlb is the squared expectation value 
of the operator at finite volume. The form factors used in (IVLlb are modified by the set of indices ui,..., (T„. The meaning of 
this index is that if (jj = — 1, the j-th incoming particle (antiparticle) is crossed into an outgoing antiparticle (particle). All the 
excitations with aj = 1 are in the incoming state. 

We would like to point out that the validity of the LM two-point function has been questioned in References lUt] and i^. 
The main concern in Ref. lIK] is that the form factors used in the formula are those found at infinite volume, and they are not 
appropriate to find finite-volume correlation functions. The problem with using infinite-volume form factors is that the energies 
are dressed at finite volume. Thus when calculating finite-volume form factors one should use the appropriately dressed form 
factors. However, as we discussed before, at large N, the TEA pseudo energies of the PCSM are trivial. The pseudo energies 
from the TEA are those of a free theory, and “undressed”, infinite-volume form factors seem appropriate. In this sense, our case 
is similar to free theories, where the LM formula is valid |@]. A similar case is that of the thermal deformation of the Ising 
model. This is a theory of free massive fermions, and the two point functions were calculated in iIitIi . 

A different objection to the LM two-point formula is discussed in i^. The authors suggest that the series (IVLlb is not well 
defined for n > 3. The form factors with both incoming and outgoing excitations have poles at real values of the rapidities. 
Each rapidity has to be integrated over the real axis, and so, the integrals in (IVLlb are divergent, and not well defined in general. 
Several regularization schemes for dealing with these divergences have been proposed ifisll .ll^. In our model, however, we will 
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see when crossing excitations to the outgoing state, the poles are not pushed towards the real axis, and all our integrals are well 
defined. If two incoming excitations have a pole at the rapidity difference 9 = iri, and one of these excitations is crossed into the 
outgoing state, the pole is moved to 6 = ±27rf, instead of 0 = 0. This is a consequence of the fact that our poles are not periodic 
under 9^9 + 2Tri, as were the usual poles considered in 

We do not have any further proof that the two-point LM formula is valid in our case, except that the usual objections against 
it do not apply. The main point we want to make is that the thermal correlation functions are not the trivial ones of a free theory. 
Even if the LM formula is not completely accurate, it is useful enough to show that the thermal correlators at large N are not 
trivial. 

We now find the general form factors needed for (IVI.ll) . Because of the global SU(A^) x SU(fV) symmetry of the PCSM, the 
non-vanishing form factors are 

{A, 9M+M' + l,bM+M’ + l,a,M+M' + l\ ■ ■ - ^A, 9M+M'+k,bM+M'+k, aM+M'+k', P, ^M+M'+k+l, dM+M'+ k+l,bM+M'+k+l', 

. . . ; P, 9M+M'+k+k' ,aM+M'+k+k' ,bM+M'+k+k' I'&boao (0) 9i,bi,ai; 

.. .;A, 9m, bM, aM',P, 9M+i,aM+i,bM+i', ■ ■ ■',P, ^m+m' ,aM+M’ ,bM+M'), (VI.2) 

with the condition k+M' — \ = k'+M. We define permutations a,T G SM+k' that take the set of numbers A= {0,..., M, M+ 
M' + k -\-1, ■ ■ ■ ,M + M' -I- fc -I- k'} to the set of numbers B = {M + 1, ■ ■. ,M + M' + k}. With this notation we can express 
the form factor (IVI.21 i as 

j^M'+k-^ , bM+M'+k+k') , (VI.3) 

We now introduce some further notation needed to write down a neat general expression for the function 
Faribi ,..., 9M+M'+k+k')- We define A^, as the subset of A, such that a(j) G {M -|-1,..., M + M'} for j G {0,..., M}, for 
all j G A^. Similarly A^ is defined such that r(j) G {M+1,..., M+M'} for j G {0,..., M},for all j G A].. We define A^ G 
{0,...,M} such that cr(j) G {M+M'+l,... ,M+M' + k},fora\lj G Al,andAl G {M+M'+k+l,... ,M+M'+k+k'}, 
such that a{j) G {M -|-1,..., M + M'}, for all j G A^. Finally we define A‘^ G {M + M' + k + l,..., M + M' + k + k'}, such 
that a{j) G [M + M' + 1,... ,M + M' + k{, for all j G A'^^. We similarly define A^, and A{., in an analogous way. For a 
given pair of permutations a, r, we define lA as the number of elements in the set A^ plus the number of elements in the set A{. 
Similarly, we define A as the number of elements in the sets A?^ and A is the number of elements in the sets A^ and A^, 
and A is the number of elements in the sets A'^ and A^.. These numbers satisfy the condition A +A + A + A = 2(M -|- fc') . 

The general form factor, found from (III.41 i by using the S-matrix and crossing symmetry is given by 


, ^M+M'-l-fe+fe') = K, 


n ~ n 

X P[ {9j — 0a{j) + 27rf) P[ {9j — 9T{j) + Stti) 
j^Al jeAl 




X P[ {9j 00-0) 27rf) P[ {9j 9r(j) Stti) 

jeAl jeA^ 


where 


iTax = 


_/ (-47r)^+'=', cr(j) ^ T(j), for all j, 
0 , otherwise 


(VI.4) 


We now substitute the exact form factors into the LM formula (IVI.lb . After some tedious but straight forward calculation, we 
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find that the finite-volume correlation function , for the operator O = ‘h/v/lV, is 

. oo 2 Z 21 —21 ——n^ «oo /*oo 

w(^)'' = 4^EEE E / de2i+i[fMh{02i+iT 


n 

j=i 


i =0 ni =0 n 2=0 rt ^=0 

{fi{dj)fi{dj+i)exp{-[tiej + ej+i)/V + ix{kj + %+i)]}]- 


{Oj - 6»j+i)2 -p 7r2 

[/i(6'i)/-i(6'i-ti)exp{-[f(ej - ej+i)/V^ + ix{kj - fc^+i)]}]^ 


n-'+n^ 

n 

n^+n^+n^ 


j=n^+n2 + l 
21 


{Oj - Oj+iY + 47r2 

[/-i(^j)/-i(^j-n)exp{[f(ej- -h ej-+i)/y -h ix(% -h %+i)]}]^ 

{Oj - 6lj+i)2 -P7r2 

[/-i(^'f)/i(^'f-ti)exp{-[f(-ej -h ej+i)/V + ix{-kj + %+i)]}]- 


n 

i+n* 

21 

n (0 0 )2_^47r2 

n2+n4 + l ^ ^ 

^ oo 21 2l—n^2l — n^—n^ «oo poo 

2 -EE E E / de,... d 02 i+Af-iidi)f-i{d 2 i+ir 


4:TT 

xn 

f=i 

X n 


l—O n^=0 n2=0 n3=0 


[/-i(^j)/i(^j+i)exp{-[f(-ej- -h ej+i)/V + ix{-kj + kj+i)]} 

{0j - 6»j+i)2 -p 47r2 

[/i(6'j)/i(^'j+i)exp{-[f(ej -h ej+i)/V + ix{kj + fc^+i)]}]- 


n^-t-n^+ra^ 


n 


21 


n 


(0j - 0j+if + 7r2 

[/i(^j)/-i(^j-n)exp{-[f(£j- - ej+i)/V + ix{kj - fc^+i)]}]^ 

(0^.-0^.^i)2+4^2 

[/-i(^'i)/-i(^'i-ti)exp{[f(ej -h e,+i)/V + ix{kj + fc^+i)]}]' 


j=n^+n2+n^ + l 


(0j - 0j+if + 7r2 


O 


(VI.5) 


We now want to study how this correlation function diverges at a; = 0. As we did in the previous section, we will take a; = 0 
and introduce a rapidity cutoff to regularize any divergence. The function (IVI.5b becomes 


w"(or = ^EE E E I"d0,... rd02i+i[h{0i)fi{02i+i)]^ 


cxD 21 2l — n^2l — n^—n‘^ 


47r 


l—O ni =0 n^—0 n'^—0 


'-A 


n 




i n^+n^ 


(0, - 0,+i)2 + 7r2 (0^ - 0^+,y + 


n 


[/i(0,)/i(0,+i)]^ 


• 1 n^+n^+n^ 


oo 21 2l — n^2l—n^—n^ pX 


■ 47r2 


n 


[/i(0,)/i(0,+i)]^ 


21 

n 


(^a - id, - 0,^,r + 4-^ 


+^EE E E I d0i... I d02i+i[f-ii0i)f-i{02i+i)]^ 


l—O n^—O n^—0 n'^—0 


f-X 


X 

-X 


n 


[f.,{0,)hi0,+,)]^ "Vt' [/i(^a)/i(^a+i)] 


i n^+n^+n^ 


n 


( 0 , - %+ l )^ + 4^2 - 0 ,+,Y + .2 ^^^11 ( 0 ^. _ 0^.^,)2 + 


n 


[/l(0,)/l(0,+l)]^ 


j—'nX +n^ + l 


47r2 


n 


j—n^ +n2+n^ + l 


{0, - 6»j + i)2 -P 7r2 


(VI.6) 


In the following section we evaluate the expression (IVI.6I) taking two very different limits. First we examine (IVI.6b at very 
large volume, V —?> oo, and recover the previous results from Section III. We then examine the opposite limit of very small 
volume, y —>■ 0. 
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VII. THE P-REGIME VS. THE e-REGIME 


In this section we examine the function (IVI. 6 I) for both very large and very small volumes. There is only one length scale in the 
PCSM, namely m. By large volume, it is meant, that V >> 1/m. In the finite-voIume-QCD literature ifl^ . this is commonly 
called the p-regime. By small volume, it is meant that V << 1/m, which is commonly called the e-regime. If instead the 
direction is made finite, the p-regime and the e-regime correspond to the low-temperature, and the high-temperature limit, 
respectively. 

The volume dependence of the expression (IVI. 6 I) is included only in the functions 


M0) 


— Vm cosh ti ; ^ 1; 

g— V m coah 6 

Vm cosh 0 ; ^ ^ 


Taking the large-volume limit, these become 


lim fa (9) 
V—^OO 


1, cr=l, 
0, cr = -l 


Then at large volume, equation (IVI. 6 I 1 becomes (IIII.71 I and we simply recover the results we calculated at infinite volume. We 
consider this limit as a trivial consistency check of the LM formula. 

Now we examine ( IVI. 6 I ) for V << 1/m. The argument we will use is similar to the one used to find the expression ( IIII.3I ). 
We use the fact that for very small V, the function ® becomes approximately a plateau, with value 1 for —C < 9 < C, 

and 0 elsewhere, where £ = In In this limit, then 

’ m 1 / ’ 


h{9) 


i, -C<9<C, 
1 , otherwise, 


Defining the new operators 


f-m 


-C<9<C, 
0 , otherwise. 


{9\i^m 


1 

2 (6» - -f 27r2 ’ 


{9\P\9') 


1 

2(6» - 6I')2-f 87r2 ’ 


{9\i\9') 


1 

(0 - 6»')2-f 7r2 ’ 


the expression (IVI. 6 I 1 becomes 




Ttt 


d9'd9{9'\ 


'-c 


i-(fi)2y \i-if^) 


|0) 


1 


-f— / d9'd9{9'\ -^ 


1 


1 


I")+1 i 


-f>V^(0). 


1 


i-(£)2y VI 


1 


l-f2 


\ 0 ) 


(VII. 1) 


The first term in the right-hand side of (IVII.lb has no dependence on A. 

We now express the result (IVII.ll) in terms of the eigenvalues and eigenfunctions of the fractional Laplacian operator. We 
switch to new variables Uj = 9j/C. For very large C (deep in the e-regime), Eq. (IVII.ll) can be written as 


W^(0)^ 




1 

l-f2 


\u) 



du(pri{u) 


2 


+ 0 {£)< 0 , 


where the operator, T, was defined in Eq. (IIII. 6 I 1 . It is important to notice that in the e-regime, the difference 5IF^(0) = 
yV^(O)^ — yV^(O) is always negative. This means that at very small volumes (or very high temperatures) the correlation 
function becomes less and less divergent at x = 0 . 
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vm. CONCLUSIONS 

We have computed the thermal expectation value of the trace of the energy-momentum tensor of the PCSM at large N, using 
the Leclair-Mussardo formula. This value is usually easy to determine from the TEA, which in our model is that of a free theory. 
The value we obtain from the LM formula is not trivial, and does not agree with what is expected from the TEA. This implies 
that there is a nontrivial ground state energy. The discrepancy arises from the fact that the two-particle S-matrix is nontrivial if 
the particles share a color contraction. The TEA only takes into account the symmetric part of the S-matrix, which is trivial at 
large N. 

We have also calculated the two-point correlation function of the renormalized field operator in a finite volume. This calcula¬ 
tion was done using the LM formula for two-point functions. The validity of this formula has been questioned before. However, 
we argued that the usual objections do not apply in our particular case. In our case, all the integrals in the LM formula are well 
defined, and pseudo energies from the TEA are those of a free theory. For very large volumes (in the p-regime), we recover 
the standard, infinite-volume two-point function from |l3t],|@], which diverges logarithmically at short distances. For very small 
volumes (in the e-regime), we saw that this divergence gets softened as we reduce the volume size. 

The ground state energy (and some excited states) of the finite-volume PCSM for some small values of N have been calculated 
before in Ref. by solving the associated Hirota equation. The large-N limit extrapolation of these results is not yet well 
understood. These results, however, only use the symmetric part of the S-matrix, so it is doubtful that our nontrivial results 
can be reproduced by simply extrapolating to large N. The authors of Ref. ll^ propose their results can be used to study the 
alternate large-A^ limit of Ref. cn. 

One might hope to obtain a nontrivial partition function and expectation values from the TEA by working at general finite 
N, and taking the large-A^ limit only at the end of the calculation. This exercise is much harder than our case, since the Eethe 
equations are highly nontrivial. However, we argued that this process should yield the same free boson partition function, as 
long as we take the large-7V limit described in this paper, and not the one from fill] 

Our interpretation of our results is that the TEA partition function, starting from the large-A^ limit of the S-matrix, is not 
enough to describe all the thermodynamics of a matrix-valued theory. There are contributions to vacuum expectation values of 
operators that arise from the matrix structure of the fields, which contains information ignored by this partition function. 
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